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(N , 

! 1 Introduction 

(N 
CN 

\^ • The aim of this work is to study some asymptotic properties of a Heath-Jarrow-Morton 

CN \ model of the term structure of interest rates driven by an infinite dimensional Levy noise. 

In particular, denoting by u{t,x), t, x > 0, the forward rate at time t with maturity 
• t + X, we shall be concerned with the following stochastic partial differential equation 

(SPDE): 

du{t,x) = [u^{t,x) + f{t,x)]dt + {a{t,x),dYo{t)), (1.1) 

^ , where Yq is a Levy process taking values in a Hilbert space with inner product (•,•), cr 

is a deterministic volatility term, and / is such that discounted prices of zero-coupon 
bonds are local martingales. Precise assumptions will be stated below. Note that in the 
special case where Yq is a Wiener process, (jl.ip can be written in the more familiar form 



X 



CXD 

du{t,x) = [u:cit,x) + f{t,x)]dt + '^a''{t,x)dwk{t), (1.2) 

k=l 



where Wk are real independent Wiener processes and / satisfies the well-known HJM 
drift condition ([l5], [lO], [H]) 

f{t,x) = TaHt,x) / a''{t,y)dy. (1.3) 
r~, Jo 



k=l 
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Invariant measures and the asymptotic behavior of ()1.2p . in the time- homogeneous case, 
have been studied by several authors (see e.g. [23], [30], [29]), allowing also the volatility 
coefficient to depend on the forward rate itself. Indeed it is widely accepted that mean 
reversion is a characteristic property of the dynamics of interest rates, and it is supported 
by empirical findings. 

On the other hand, the literature on HJM models driven by Levy processes has 
considerably grown in the last few years: let us just cite, among others, [Ij, [8j, and 
of course the work from where it all started [2] (where general random measures are 
added to Brownian motion as driving noises). The asymptotic behavior of such models, 
however, does not seem to have been addressed. The present paper offers a first step in 
this direction, in a simple model with deterministic volatility. The setting is similar to 
that of [30j (where the case of HJM models driven by Wiener process was considered), 
but the choice of state space is different, and the well-posedness of the model is somewhat 
more complicated. 

The paper is organized as follows: in section [2] we derive sufficient conditions on / 
ensuring that the bond market is arbitrage free and write the SPDE (jl.ip as an abstract 
evolution equation in a suitable Hilbert space of forward curves, about which we discuss 
existence and uniqueness of mild and weak solutions. In section [3] we discuss existence 
and uniqueness of invariant measures, as well as the convergence in law of forward curves 
as time goes to infinity. 

After the first draft of this paper was completed, the author was informed of some 
papers (most of which not yet published nor posted to any standard preprint server) 
with some overlap with the present one. In particular, drift conditions were derived in 
[24| and in [T7]. Our derivation, as well as the setup, is slightly different, and we include 
it for the sake of completeness. Existence and uniqueness of local mild solutions for HJM 
models driven by finitely many independent Levy processes is discussed in [11] (see also 
[21]), where the volatility is also allowed to be state dependent. In the present paper 
we consider a possibly infinite dimensional driving Levy process, which is not supposed 
to be the superposition of independent one-dimensional Levy processes, and we look 
for global solutions. Moreover, in [5S] the authors obtain conditions for existence and 
uniqueness of global mild solutions with an infinite dimensional Levy process and state 
dependent volatility. However, their state space is a weighted L2 space, which seems 
inappropriate for modeling purposes: in fact very irregular forward curves cannot be 
excluded. We would also like to point out that the choice of state space in |11] is the 
same as in our paper, and for this reason the former paper only treats local solutions. 
On the other hand, as already observed (see also [lOj for an extensive discussion), our 
choice of state space is better than that in [25]. Finally, the existence and uniqueness 
statement of proposition [8] can be deduced by the result in [27] (see also [28] for a result 
on ergodicity of HJM models with one-dimensional, possiblly multiplicative noise, in 
spaces of weighted L2 curves). 

Let us conclude introducing some notation. Given two separable Hilbert spaces H, 
K we shall denote by C{H,K), Ci{H,K) and C2{H,K) the space of bounded linear, 
trace-class, and Hilbert-Schmidt operators, respectively, from H to K. Cf stands for 
the subset of £1 of positive operators. We shall write Ci{H) in place of Ci{H,H), 
and similarly for the other spaces. Given a self-adjoint operator Q € C^{H), we set 
\x\q := {Qx, x), X H. The Hilbert-Schmidt norm is denoted by | • I2. The characteristic 
function of a set A is denoted by XAi and Xr stands for the characteristic function of 
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the set Br := {x £ H : \x\ < r}, where H is a Hilbert space. Given a continuously 
differentiable increasing function a : M-|- [^,00) such that a~^^^ € Li(M-i-), we define 
^2 a '■— ^2aO^+) the space of distributions (f> on M_|_ such that \(j)^"\x)\'^ a{x) dx < 
00. 



2 Musiela's SPDE with Levy noise 

Throughout the paper Yo{t), t > 0, shall denote a Levy process taking values in a (fixed) 
Hilbert space K, with generating triplet [bo, Ro,mo], i.e. 

logEe^<^'^°«) =i(6o,y) - l{Roy,y) + / (e*^«'^> - 1 -^(e,y)xi(0)"^o(dO, 

with bo £ K, Rq £ Cf{K), and mo a a-finite measure on B{K), the Borel cr-algebra of 
K, satisfying 

"io({0})=0, /" |e|'mo(dO<oo 

(see e.g. flS] for details). This integrability assumption serves two purposes: it ensures 
that E|yo(t)P is finite for all t > 0, thus allowing to construct mild solutions of the 
SPDE (jl.ip via an L2 theory of stochastic integration, and it allows to use Fubini's 
theorem to establish no-arbitrage sufficient conditions. The assumption will turn out 
not to be a real restriction, as the no-arbitrage condition essentially requires existence 
of exponential moments of mo (see also [Ij). 



2.1 Drift condition 

Let us consider, in the spirit of the original paper [15] (see also [14]), the following 
integral equation 

u{t,x) = u{<d,x) + [ b{s,x)ds+ [ {a{s,x),dYo{s)), (2.1) 
Jo Jo 

where b and a are random vector fields predictable in s and Borel measurable in x. In 
particular, they could depend on u itself. 

We shall give conditions under which the dynamics (j2.ip is compatible with a no- 
arbitrage hypothesis, namely that the corresponding discounted bond prices are local 
martingales. The arguments of the proof consist of the Levy-Ito decomposition in Hilbert 
space (see e.g. ft]) and the calculus for random measures, following j^. We denote by 
P{t,T), <t <T, the discounted price of a zero-coupon bond expiring at time r. 

Theorem 1. Let T,(t,x) = — a{t,y) dy and E{t,x,^) = — {a{t,y),£) dy for all 
x > 0, £ K . Assume that for all t < 00 

rt roo ft roo 

/ / \b{s , x)\ dx ds < 00 , / / \a{s,x)\'^ dx ds < 00. (2-2) 
Jo Jo Jo Jo 

Moreover, assume that for all x > one has 

1 

T 



b{t,y)dy = u{t,x)-u{t,0) + -|S(t,x)||^ 



+ f (e^(^^-^^^-l-E{t,x,Oxi{0)mo{dO 

J H 



(2.3) 
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dF X dt-a.e. Then the discounted bond price process t ^ P{t,T), t > 0, is a local 
martingale for all t >t. 



Proof. The Levy process Yq admits the decomposition 

/ ^N{ds,dO+ f 



Yo{t) = bot + W{t) + f I i N{ds, dO + f I e N{ds, dO, (2.4) 
Jo J\£\<1 Jo J\£\>1 



where bo (z K, W is a valued Wiener process with covariance operator Rq, is a 
Poisson measure on K with compensator mo, and N(ds, d£,) := N(ds, d^) — (ismo(dO- 
Using the decomposition (j2.4p one can write 



u{t,x) = u{0,x)+ / b{s,x)ds+ / {a{s,x),dW{s)) 
Jo Jo 

+ / {cr{s,x),dz{s)) + / {a{s,x),dz{s)), 



where ^ ^ 

m= [ I iN{ds,di), z{t)= [ [ CN{ds,d^). 
Jo J\£\<i Jo J\£\>1 

However, let (ek)keN be a base of K and set x'^ = {x, e^), x (z K. Then one has 

t-t " rt 

/ {a{s,x),dz{s)) = lim / a''{s,x) dz''{s) 

Jo Jo 

= lim V /" /" cj''(s,x)e''iV(ds,de) 

= hm / / j;a'=(s,x)e'^iV(d.,de), (2.5) 



where in the second line we have used the associativity of stochastic integrals (i.e., 
using the notation of the stochastic calculus of semimartingales and stochastic measures, 
H-{K*n) = {HK)*fi - see e.g. [E]). Since Ylk=i o'^(s,x)^'^ converges (for any fixed x) 
to (cr(s, x), in L2{^ x [0, t] x H, V, d¥ x ds x drrio) for all t > 0, with V the predictable 
c7-field, one finally obtains 

{a{s,x),dz{s)) = / rj{s,x,C)N{ds,dS,), 



Jo J\£\<1 

where rj(s,x,(^) = {a{s,x),S^). Similar reasoning shows that the same type of identity 
holds for integrals with respect to A^. 

For X > 0, let us denote by p{t, x) the discounted prices at time i of a risk free zero- 
coupon bond expiring at time t + x. By the definition of p{t, x) and the equation for u 
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we obtain 

/•t PX 

logp{t,x) = - u{s,0)ds- / u{t,y)dy 
Jo Jo 

ft fX 

u{s,0)ds— / u{0,y)dy 
lo Jo 

b{s,y)dsdy- / / {c7{s,y),dW{s))dy 
Jo Jo 

-ril v{s,y,ON{ds,dOdy- r [ [ rj{s,y,ON{ds,dOdy 
Jo Jo J\^\<i Jo Jo J\^\>i 

ft fX ft ft 

- / ^(5,0)^5- / u{0,y)dy+ / B{s,x)ds+ / {T.{s,x),dW{s)) 
Jo Jo Jo Jo 



+ / E{s,x,ONids,dO+ / E{s,x,ONids,dO, 
Jo J\(,\<i Jo J\(,\>i 

where the third equahty follows by Fubini's theorem (see [5] and [18j ) and the definitions 
of S and E, together with B{t, x) := — Jq b{t, y) dy. 

Applying Ito's formula, setting Ci^^x) = logp{t,x), one gets 

p{t, x) = e'^^*'^) = e^(°'^) + e'^(^-'^) ( - u{s, 0) ds + B{s, x) ds + (S(s, x),dW{s))'^ 



+ 1 e<^{s-,x) ^[^^ QC^^^ X) + f f (^e<is-,x)+Eis,x,0 _ ^({s-,x)^ ^(^^^ 

2 JO Jo J\i\<i 

+ f [ (e^(^-'^)+^(^'^'«) - e^(^-'^)) N{ds, dO 
Jo J\^\>i 

Jo J\i\<i 



or equivalently 



dp{t,x) _ ( ^ ^ 1,^^^ ^ 2 



n(t, 0) + Bit, x) + - dt + x),dWit)) 



p{t-,x) V ^ ' ' ^ ' ' 2 



+ 


/ ^ 


'^f,E{t,x,0 


- l)N{dt,di) 




/|CI<i 






+ 


/ ^ 


'^f,E{t,x,0 


- 1) N{dt, di) 




/|CI>i 






+ 


/ ^ 


'^f,E{t,x,0 


- 1 - E(t,x,(,))mo{d£,)dt 




/lfl<l 







For T > t, by r) = p(t, r — t) it follows that dP{t, r) = dp{t, t — t) — px{t, t — t) , and 

p{t,T-t) = e-/o"(*'°)'^"P(t,T-t) implies p^(t,T-t) = e" io* r - 1). Setting 

P(i,r) = e-^o<^^y)<^y, one has Px{t,T - t) = -P{t,T - t)u{t,T - t) and Px{t,T - t) = 
—u{t, T — t)p{t, T — t), and finally 

dP{t, t) = dp{t, T-t) + u{t, T - t)P{t, t) dt. 
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Together with the equation for p{t,x), this imphes 
dP{t, t) 



-u{t,0)+uit,T-t) + B{t,T-t) + l\j:{t,T-t)\l^dt 
P{t-,T) V ^ ^ 

+ {^{t,T-t),dW{t))+ [ (e^(*'^-*'«) - 1) iV(dt, de) 
+ / (e^(*'"-*'«) - 1) N{dt, dO 
+ / (e^(*'-*'f) - 1 - E{t, T-t, 0) moidO dt. 



'l?l< 

Lightening notation a bit, one can write 



/ {e^-l)dN= [ {e^-l)dN+ [ {e^ - 1) dmo dt , 
J\£\>1 J\£\>1 J\£\>1 



I5I>1 ■''ICI>1 Jm>^ 

hence P{t, r) is a local martingale if 



= -u{t,0)+u{t,T-t) + B{t,T-t) + -\J:{t,T-t)\l^ 



+ I (e^(*'^-*'«)-l)moK)+ / (e^(*'^~*'«)-l-i^(t,r-t,0)mo(de) 
J\i\>i J\i\<i 



= -u{t, 0) + n(t, T-t) + B{t, T-t) + r - i) 

Jk 

and the theorem is proved. □ 

Remark 2. The above theorem implies a "drift condition" that generalizes the HJM 
condition (jl.3p . In particular, assume that Y{1) admits exponential moments, or equiv- 
alently that J^j e^^'^^ mQ{d(^) < oo for all z & H, and define the function 'iIj{z) = 
log Ee^^'^(^)^ Then taking into account and (|2.1|) . (j2.3p implies the following 

relation between the drift and the volatility functions: 



f{t,y)dy = ip(- / a{t,y)dy) (2.6) 
^ Jo ' 



for all t E [0,r]. Unfortunately this identity is "implicit", and only under further 
assumptions can it be made more explicit (see below). 

2.2 Abstract setting and well-posedness 

We shall rewrite the SPDE (jl.ip as an abstract stochastic differential equation in the 
space H = L\^. The space H endowed with the inner product 

{(p^il)) = I (p' (x)ip' (x) a{x) dx + lim (p{x)ip{x) 
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is a separable Hilbert space. This choice of state space is standard and is apparently due 
to Filipovic [To]. Nonetheless, other authors have studied related SPDEs in different 
function spaces, e.g. in weighted L2 spaces, weighted Sobolev spaces, or fractional 
Sobolev spaces (see [H], [30], pj respectively). 

Let us define on H the operator A : / 1— > /', with domain D{A) = L\ L\ ^, which 
generates the semigroup of right shifts [e*^(/)](x) := ^(x + t), t > 0. Musiela's SPDE 
(jl.ip can be written in abstract form as 

du{t) = {Au{t) + f{t)) dt + B{t) dYo{t), (2.7) 

where f{t) = f{t,-) and B{t) G C{K,H) is defined by [B{t)u]{-) = {a{t,-),u)K, with 
suitable regularity assumptions on a. 

Several papers deal with the solution of this type of equations in the time-independent 
case with f{t) = 0, B{t) = B. Here we limit ourselves to mention [3j, which is proba- 
bly the first paper considering weak solutions (in the sense of PDEs), and [12], where 
an analytic approach is used to solve, even in the strong sense, equations of the type 
(j3.ip . possibly in a larger space, allowing the characteristic function of Yq{1) to be only 
continuous instead of Sazonov continuous (see also [19\, ^20j, [26j). 

Our goal in this subsection is less ambitious, namely we shall only prove that the 
formal solution of (12.70 given by the variation of constants formula 

n(t) = e*^uo+ f e^'~'^^f{s)ds+ f e^'~''>^B{s) dYo{s) (2.8) 
Jo Jo 

is a well defined process and provides the unique weak solution to (|2.7p . In particular, 
this implies that the Levy-based model (II. ip for the evolution of forward curves is well 
posed under appropriate assumptions on a. 

In analogy to [i], we shall prove that ()2.8p belongs to the space 'H2{T) of mean 
square continuous process on [0, T] with values in H adapted to the filtration generated 
by Y . We endow 7^2 (^) with the norm defined by 

||F||2 := sup nF{t)\]j. 

te[o,T] 

Proposition 3. Assume that 



[ {\f{t)\H + \B{t)\l)dt<cx^. 
Jo 



Then u defined as in f2. 8\) belongs to T-L2{T). 

Proof. Adaptedness is immediate by definition. Using the Levy-Ito decomposition in 
the form 

t 



Yo{t) = at + W{t)+ [ [ CN{ds,dO 
Jo Jk 



with a = 60 + i|^|>i C "^o('^O) ^"^^ taking into account Propositions 2.2, 2.3 in it is 
enough to consider the case when YQ{t) has no drift and no Brownian component. In 
particular one has 



E 



f e^t-')AB{s) f ^N{ds,di) <N^ [ \B{t)\l [ {^{"^ mo{dC) dt < 00 
Jk Jo Jk 
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with = suptg[o^T] k*^!- Moreover \ J^^ e^^~^^^f{s) ds\ < N jQ\f{s)\ds < oo, hence 

||u||2 < OO. 

Let us now prove that tH^E|ti(t) pis continuous. Setting YA(i) := e^^'^'I'^Bir) dYo{r), 
it is enough to prove that t i-^ E|y4(t)p is continuous. For < s < t < T, write 

YA{t)-YA{s)= r {e^'-'^-^^-e^'-'^^)B{r)dYo{r)+ f e^'-'^^B{r) dYo{r), 





where the two terms on the right-hand side are uncorrelated. Since 

2 



E 







(e(*-0^-e(^-'^)^)S(r)dyo(r) 



< r \e^*-^^^ - e^'-^^^\^\B{r)\ldr I |eP"io(dO^O, 



K 



E 



'~i-^)^B{r)dYQ{r] 



< 



* \Jt-r)A\2 



B{r)\ldr [ |e|2mo(de) ^0 
JK 



as s t i, the result fohows. The case < t < s < T and s j t is completely analogous, 
hence omitted. □ 

The following proposition shows that existence in the mild sense for equation ()2.7p 
implies existence and uniqueness in the weak sense. This fact was essentially proved by 
A. Chojnowska-Michalik in ^ in the time-independent case. Here we give a more direct 
proof that closely follows [5j . 

Proposition 4. Equation (^7}) has a unique weak solution given by 112. <^) . 

Proof. Let us define the additive process Y{t) = f{t) + £B{s)dYo{s) and write (pTll 
in the form 

du{t) = Au{t) dt + dY{t), u{0) = uq. 

It is enough to consider the case uq = 0, the extension to the general case being imme- 
diate. We need to prove that, for v € D{A*), 

\^'-''>^dY{s),v) = f { r e^'-''''>^dY{r),A*v)ds + {Y{t),v), 
^0 Jo 

or equivalently 

JO Jo Jo 

In fact, the right-hand side of the previous expression can be equivalently written as 

A*v,dY{r))ds, 



t rs 



Jo 



hence, using the identity 

ft-S 



qTA* a* 







A*vdr = (e(*-")^* - l)v, 
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Pubini's theorem yields the desider conclusion. Uniqueness will follow if we prove that 
a weak solution is a mild solution. One immediately recognize that the proof of Lemma 
5.5 in [5], repeated word by word, yields the identity 

{w{t), cp{t)) = f\w{s),^'{s) + ds + f\<P{s), dY{s)), t G [0, T], 



where w is a weak solution of and G C^([0, T],D{A*)). Taking 0(s) = e(*"")^*0o, 
00 G implies 

{w{t),M = { fe^'-'^^dY{s),4>o), 
Jo 

hence u{t) = w^t) because D{A*) C H densely. □ 

It is clear that in order to obtain results on the well-posedness of the Musiela's 
SPDE (II. ip it is necessary to establish conditions on a and Yq such that the hypotheses 
of Proposition [3] are satisfied. The following sufficient conditions are obviously the most 
general ones, but also the hardest to verify. 

Proposition 5. Let {ek)keN be a basis of K, and define a^{t,x) = {a{t, x), ek)K ■ As- 
sume that 

(■T foo 

(2.9) 



pi ^ poo 

/ VJ / {t, x)\'^a{x) dx dt < oo 
Jo fc^^Jo 



and that there exists / : [0, T] ^ satisfying h2. 6)) such that \f{t)\ dt < oo. Then 
there exists the mild solution of [27. 



Proof One has \B{t)\l = J2T=i \B{t)ek\'^, and B{t)ek = a^{t, •), hence 



r-oo 

\B{t)ek\^ = I \a^{t,x)\^a{x)dx. 
Jo 



Therefore ^Ml implies /g \B{t)\ldt < oo, and the result follows by proposition [3l □ 

Much simpler conditions can be stated if the driving noise Yq is finite dimensional. 

Proposition 6. Assume that Yq is a W^-valued Levy process such that Jjg„ e^^'^^ rno{d^) < 
oo for all z G Br, for some r > 0. // Jq a{t, y) dy G Br for all {t, x) G [0, T] x M+, and 

pT poo 2 

J [J \a^{t,x)\^a{x)dxj dt < oo (2.10) 
for all k = 1, . . . ,n, then \2. 6)) reduces to 

f{t,x) = -(^a{t,x),Di;(^- r a{t,y)dy^y (t, x) G [0, T] x M+, (2.11) 



and {2. T\ j admits a mild solution. 
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Proof. It is enough to check that the hypotheses of proposition [3] are satisfied. In 
particular, similarly as before, one has, using Jensen's inequality, 



\B{t)\ldt= y2 \cr^{t,x)\'^a{x)dxdt 
Jo ^^^Jo 



fc=i"'" k=i 



where the last inequality is immediate by (j2.1Up . Moreover, in (|2.11|) the quantity Dip{x) 
is well defined for |x| < r because ip € C°°{Br). This implies that there exists a positive 
constant N such that < N, \ip^i^3{z)\ < N for all z G Bj-- We have 



fx{t,x) = -(^aa:{t,x),Dip(^- J a{t,y)dyYj 

+ / cj{t,y)dyy{t,x),a{t,x)), 



hence \f{t)\Hdt < oo \i J^[a'^{t, x)a^ {t, x)]'^a{x) dx dt < oo for all i, j < n. The 
condition a{t, y) dy < r for all x > implies that limy^oo Wit, y)| = for all t G [0, T], 
hence, using Cauchy-Schwarz' inequality, 

a^{t, x)a^ {t, x)]'^a{x) dx dt 

T foo V 1 /9 / rT roo 



JQ 



~U J '^*(*'^)^"(^)^^^*)^^^(_/| (J^{t,x)^a{x)dxdty^^ 
<iV(a)( r \a\t,-tHdt)"\ r\a\t,-%dtf^ 



where the second inequality follows by (5.8) in [TO]. □ 

An analogous expression could be obtained for a general Hilbert space valued noise 
Y , if one can guarantee that ip is Frechet differentiable. In the next proposition we shall 
identify the Frechet derivative D'il){x) G >C(i^, R) with its Riesz representative vector in 
K. 

Proposition 7. Assume that e^^'^^ mo((i^) < oo for all z € Br, -0 ^ C^i^r), for 
some r > 0, and 

/ (y]/ Wlit, x)\'^ a{x) dx] dt <oo (2.12) 
Jo \_,7o ^ 

If Jo o"(*i y) dy G Br for all {t, x) G [0, T] x R+, then fO) reduces to ^KTT]) . where (•, •) 
is the inner product of K, and equation (2.1) admits a mild solution. 

Proof. As seen before, we have 



B{t)\l = Y. / \al{t,x)\^a{x)dx, 
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hence, by (f2J2]) . \B{t)\l dt < (T \B{t)j dtf^"^ < oo. Let us prove that \fit)\dt < 

00. The same expression for fx{t,x) as in the proof of the previous proposition holds, 
mutatis mutandis. Therefore, in view of (f2T2]) and ^p G C|(S^), it is enough to show 
that Jq \a{t,x)\^a{x) dx dt < oo. Let (cfc) be a basis of K, and set an{t,x) = 

Y^=i ^^{ii ^)^k- Let (j)n^ a smooth approximation of x i— > |x| in M" such that \D(t)i\ < 

1, then we have as follows by (5.8) in [TU] 

) 

(t)^^\an{t,x)fa{x)dx < Ar|0(f)(a„(t,.))|J, 

= N[j^ \D,,<p^^\<jn{t,x))\^a{x)dx) 

\ 2 

\Dx(Jn{t,x)\^a{x) dxj 

oo 2 

<^(E / W'it,x)\'aix)dx) 
k=i 

with N = iV(a), thus by (l2T^ 

/ / (t)^a\'^n{t,x)fa{x)dxdt<Nl (V/ |cT^(t,x)pa(x)dx) < oo. (2.13) 
./o JO Jo ^ i^^^JO ' 

Since the bound in (j2.13p does not depend on e nor on n, passing to the limit as 
e — > we get Jq |o"n(i, x)|'*a(x) dx < oo, and letting n tend to infinity we finally get 
So Io° |o"(t, x)l^a(x) dx < oo. □ 



3 Invariant measures and asymptotic behavior 

In this section we assume cr(t,-) = (t(-), thus also B{t) = B, f{t) = /. In view of 
the no-arbitrage considerations in the previous section, we also assume that mo admits 
exponential moments, and B G C{K, Hq), where Hq := {g ^ H : g{oo) = 0}, hence also 
/ G Hq. In order for the following results to hold, it is not necessary to assume that / 
is such that no-arbitrage is verified, even though this is of course the situation we are 
interested in. 

Let us rewrite (12. 7|) . for convenience of notation, in the more compact form 

du{t) = Au{t)dt + dY{t), (3.1) 

where Y{t) := ft + BYQ{t). Then one can easily prove that y is a i^-valued Levy process 
with triplet [b, R, m] , where 

b = f + Bbo+ [ B^{xi{BO - XiiO) moidO 

JK 

R = BRqB* 
m{di) = mo{B-^di) 

The following proposition gives a simple sufficient condition for the existence and 
uniqueness of an invariant measure for an HJM model with deterministic volatility and 
Hilbert space valued Levy noise. The only real requirement is that the state space L\ ^ 
is chosen with an exponentially growing weight a. 
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Proposition 8. Assume that a{x) = e"^, a > 0, and the forward curve at time zero is 
deterministic. Then there exists a unique invariant measure for 13. 1\) to which the law 
of u{t) weakly converges as t ^ oo. 

Proof. Writing equation ()3.ip in mild form, recalling that the range of B is contained 
in Hq, one recognizes that u{t, oo) = uo(oo) for all t > ("long rates never fall"). 
Considering the isomorphism H = Hq © M, ()3.ip is equivalent to the system 

f du{t) = Au{t) dt + dY{t) 2^ 
1 oo) = ^, 

where u is the projection of u on i^o, A still denotes the restriction of A to Hq, and 
^ G R. Let us show that e*^ is exponentially stable on Hq: 

\e^^(l)\%^ = / (P'{x + tfa{x) dx < e~"* / (l)'{xfa{x) dx = e-"*|(/.||^g, 

JO 

i.e. le*"^! < e~*°/^. The obvious inequality > log(l + x), x > 1, and the assumption 
Ik\^\^ mo((i^) < oo imply that 



/ log(l + lel) m((iO < / \i?mQ{B-^di) 

< |Sp / |CP"ioK) < oo. 
Jk 



Therefore theorem 6.7 of [H] yields the existence of an invariant measure p, on Hq for the 
first equation of ()3.2p . hence fi = p, ^ 6e is an invariant measure for (jS.ip on H. Since 
|e*^| < e~*°/'^ — > as t ^ oo (i.e. e*"^ is stable), proposition 6.1 of [3j (or theorem 3.1 
of [l2] ) imply that fl is infinitely divisible with triplet [boo , Roo ; fnoo] , where 



6oo = lim / e'^bds + / / e'^^ixiie'^0 " XiiO) MdO ds (3.3) 

L 7o Jo J Ho J 

/■oo 

i?^= / e'^^Re'^' ds (3.4) 

mM)= m{{e'^)-^di)ds, moo({0}) = 0, (3.5) 
Jo 

are all well defined thanks to the stability properties of e*^. In particular Ji is unique. 
Finally, by proposition 6.1 of [3] we have that Jl coincides with the law of the random 
variable e'^^ dY{s), and lemma 3.1 of [3] allows to conclude that the law of u{t) 
weakly converges io Jl. □ 

Remark 9. The decomposition H = Hq ©M was already used in ^29j, but essentially 
the same "trick" already appeared, perhaps less explictly, in [HU]. The fast growth 
at infinity of the weight a was needed in [29] as well, where it is assumed that a := 
\n.ix>Q Oi' {x) / a{x) > 0. In fact, Gronwall's lemma immediately yields that this condition 
implies a{x) > a(0)e"^'. 
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The choice of the weight function a, as just seen, determines the stabihty properties 
of the semigroup e*"^. For a generic choice of a we cannot guarantee exponential stabihty 
of e*"^ in Hq, but we still have stability, in the sense that \e^^g\Ho — > as t ^ oo for 
any g € Hq. However, in order to obtain existence of an invariant measure for ()3.1I) . 
the conditions to verify become quite difficult, in general. In particular the following 
characterization holds, the proof of which follows [3] or |12j . 

Proposition 10. Assume that the forward curve at time zero is deterministic. The 
following conditions are sufficient and necessary for the existence of a (unique) invariant 
measure fi for 113.1]) : 



(i) supTr / e'^BRoB*e'^' ds; 

i>0 JO 



(ii) / / (|e*^xp A 1) m{dx) dt < oo; 

Jo J Ho 

fiii) the limit in 13. 3\) exists. 



Moreover, ^ = ^jl ® 5i, where /i is infinitely divisible with triplet \boo^Roo,fnoo\ given by 
\3. 3(j - [37^) . Finally, the law of u{t) weakly converges to p, at t ^ oo. 



Proof. The semigroup e is stable on Hq because 

g\xfa{x-t)dx< g'{xfa{x)dx^^O, (3.6) 



\e 

where the inequality follows by monotonicity of a and the limit is zero because the 
integrand is in Li(M+). Therefore theorem 6.4 of [3j (or theorem 3.1 of [12j) implies 
that the infinitely divisible measure p, on Hq with triplet [boo , Roo , 'm'oo] is invariant 
for the first equation of ()3.2p . The proof is then completed exactly as in the previous 
proposition. □ 

Corollary 11. Assume that the forward curve at time zero is deterministic and that 
there exists a function (f) G Li(M_|_) nL2(lR+) such that |e*^x|//(, < 0(t)|x|/f(,. Then fi 
as defined in the previous proposition is the unique invariant measure of 113. 1\) and it is 
ergodic. 

Proof. Since e*"^ is a stable semigroup on Hq, it is enough to verify hypotheses (i)-(iii) 
of the last proposition. We have 

oo /> /"OO 

tAi2 



(|e*^x|^ A l)m(dx)dt < / / (t){ty\x\^ m{dx) dt 
Ho Jo Jho 



poo p 

< / (t){tfdt\B\^ I xp mQ{dx) < oo, 
Jo Jk 



because |xp mo((ix) < oo. Since tuq admits exponential moments, then \x\ mQ{dx) < 
oo, and 



/ / |e*^x(xi(e*^x) -xi(x))|m((ix)dt < 2 / (t){t)dt\B\ \ x\mQ{dx) < oo. 
10 Jho Jo Jk 
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Let us now prove that linit^oo Jq e^^bds exists in Hq: we have 



thus b'{x) 



b{x) := 
-b{x) and 



hm 

t— >oo 



e'Hds 



(x) 



%ds < \b\H, 





i.e. 600 £ Hq. Similarly we have 



b{s) ds, 



(s) ds < 00, 



a; > 0, 



sup 

t>o Jo 



e*^Si?oS*e^^* ds < \B\^TiRo 



(s)'^ ds < 00, 



i.e. Rno is well defined. 



□ 



Remark 12. The semigroup e*"^ is a contraction semigroup for any choice of a, e.g. by 
(13. 6p . hence if (p as in the previous corollary exists, then one can always choose \(j){t) \ < 1, 
thus (p € Li(R4_) also implies (j) ^ -^2(1^+)- A possible choice of (j) (although very rough) 
could be 

/ a(x) \i/2 
cPit) = sup , 
x.>o \ a{x + t)/ 

provided that sup^>o [a{x)/a{x + t))^^"^ dt < 00. 

We have seen in the previous proposition that one of the necessary and sufficient 
conditions for the existence of an invariant measure is that b^o exists. In particular, if 
mo is symmetric, the problem reduces to proving that b := e'^^bds is a well-defined 
element of Hq. In fact, if mo is symmetric then m is symmetric as well, and the second 
summand on the right hand side of ()3.3p is zero for all t > 0. It is thus natural to look 
for conditions on a such that the norm of b can be bounded in terms of the norm of 
b. This is indeed possible, and one can give a sharp condition, namely ()3.7p below is 
necessary and sufficient for \b\ < N\b\ to hold. 



Proposition 13. Assume that 



sup 
x>0 Jo 



a{y) dy 



dy < 00. 



(3.7) 



Then e'^^bds exists in Hq. 



Proof. As in the proof of corollary II H we only need to prove that b'^a G Li(M+). Let v 
be a nonnegative Borel measure on By a result of Muckenhoupt [22], we have that 
the following weighted Hardy inequality holds for all measurable functions / 



f{y)dy u{dx)<N / f{xfv{dx), 



where A'^ is a positive constant that depends only on a, if and only if 



sup[z^([r, +00))] 

r>0 



1/2 



^JO 



du\-^ 



dx 



dx 



1/2 



< 00. 



14 



By the change of variable x i— > x ^, choosing v{dx) = a{x) dx, we obtain that ()3.7p is 
necessary and sufficient for 



b'{y) dy a{x) dx < N b'{x)'^a{x) dx. 



Since b E Hq, we have that \J^b'{y)dy\ = \b{x)\ and that the right-hand side of the 
previous inequahty is finite. □ 



4 Conclusions 

We have considered an equation of HJM type driven by a Levy process taking values in 
a Hilbert space, obtaining sufficient conditions for the absence of arbitrage. Assuming 
that the volatility operator is deterministic, and using Musiela's parametrization, one 
obtains a stochastic evolution equation of Ornstein-Uhlenbeck type. We have discussed 
existence of mild solutions and existence, uniqueness, and ergodicity of invariant mea- 
sures, generalizing previous work of Vargiolu [HDj, who considered the situation of a 
driving Brownian motion and used a different state space. The choice of the state space 
L2 ^ seems to be the standard by now, since its elements enjoy most desirable features 
for a forward curve. If the weight function a grows exponentially at infinity, the HJM 
dynamics admits a unique invariant measure. A similar conclusion was obtained by 
Tehranchi in [29j, where a Musiela equation with state-dependent volatility and Brow- 
nian noise was considered, and the results of [6j could be applied. It would be natural 
to consider also in the setting of Levy noise a state-dependent volatility operator, but 
unfortunately there are comparatively very few results on evolution equations with jump 
noise, which need to be established first. 
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